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ABSTRACT 

We present analysis of the shape and radial mass distribution of the E4 galaxy 
NGC 3923 using archival X-ray data from the ROSAT PSPC and HRI. The X-ray 
isophotes are significantly elongated with ellipticity e x = 0.15(0.09 — 0.21) (90% con- 
fidence) for semi-major axis a ~ 10hy kpc and have position angles aligned with the 
optical isophotes within the estimated uncertainties. Applying the Geometric Test for 
dark matter, which is independent of the gas temperature profile, we find that the 
ellipticities of the PSPC isophotes exceed those predicted if M cx L at a marginal sig- 
nificance level of 85%(80%) for oblate (prolate) symmetry. Detailed hydrostatic models 
of an isothermal gas yield ellipticities for the gravitating matter, e mass = 0.35 — 0.66 
(90% confidence) , which exceed the intensity weighted ellipticity of the i?-band optical 
light, (e R ) = 0.30 (e% ax = 0.39). 



We conclude that mass density profiles with p 



are favored over steeper 



profiles if the gas is essentially isothermal (which is suggested by the PSPC spectrum) 
and the surface brightness in the central regions (r Ss 15") is not modified substantially 
by a multi-phase cooling flow, magnetic fields, or discrete sources. We argue that 
these effects are unlikely to be important for NGC 3923. (The derived e mass range is 
very insensitive to these issues.) Our spatial analysis also indicates that the allowed 
contribution to the ROSAT emission from a population of discrete sources with oc 
Y>r is significantly less than that indicated by the hard spectral component measured 
by ASCA. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: fundamental parameters 
- galaxies: individual (NGC 3923) - galaxies: structure - X-rays: galaxies. 



1 INTRODUCTION 



The structure of the dark matter halos of galaxies provides 
important clues to their formation and dynamical evolution 
(e.g. Sackett 1996; de Zeeuw 1996, 1997). For example, in 
the Cold Dark Matter (CDM) scenario (e.g. Ostriker 1993) 
there is evidence that the density profiles of halos have a 
universal form essentially independent of the halo mass or 
firj (Navarro, Frenk, & White 1997; though see Moore et al. 
1997). The intrinsic shapes of CDM halos are oblate-triaxial 
with ellipticities similar to the optical isophotes of elliptical 
galaxies (e.g. Dubinski 1994). The global shape of a halo 
also has implications for the mass of a central black hole 
(e.g. Merritt & Quinlan 1997). 



available (e.g. Sackett 1996; Oiling & Merrifield 1997)Q Stel- 
lar dynamical analyses that have incorporated the informa- 
tion contained in high order moments of stellar velocity pro- 
files have made important progress in limiting the uncer- 
tainty in the radial distribution of gravitating mass arising 
from velocity dispersion anisotropy (Rix et al. 1997; Gerhard 
et al. 1997). However, as indicated by the paucity of such 
stellar dynamical measurements, the required observations 
to obtain precise constraints at radii larger than ~ R e are 
extensive, and the modeling techniques to recover the phase- 
space distribution function are complex. It is also unclear 
whether this method can provide interesting constraints on 
the intrinsic shapes since only weak limits on the range of 
possible shapes have been obtained from analysis of velocity 
profiles out to ~ 2 R e (e.g. Statler 1994). 



At present accurate constraints on the intrinsic shapes 
and density profiles of early-type galaxies are not widely 



* The distribution of dark matter in spiral ga laxies is also far 
from being a solved problem — see, e.g. Broeils (L997). 
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Interesting measurements of the ellipticity of the gravi- 
tating mass have been obtained for two Polar Ring galaxies 
(Sackett et al. 1994; Sackett & Pogge 1995) and from statis- 
tical averaging of known gravitational lenses (e.g. Keeton, 
Kochanek, & Falco 1997), but owing to the rarity of these 
objects it is possible that the structures of their halos are not 
representative of most early-type galaxies. Moreover, gravi- 
tational lenses, which are biased towards the most massive 
galaxies, only give relatively crude constraints on the ellip- 
ticity and radial mass distribution for any individual sys- 
tem and only on scales similar to the Einstein radius (e.g. 
Kochanek 1991). 

The X-ray emission from hot gas in isolated early-type 
galaxies (Forman, Jones, & Tucker 1985; Trinchieri, Fab- 
biano, & Canizares 1986; for a review see Sarazin 1997) prob- 
ably affords the best means for measuring the shapes and 
radial mass distributions in these systems (for a review see 
Buote & Canizares 1997b; also see Schechter 1987 and the 
original application to the analogous problem of the shapes 
of galaxy clusters by Binney & Strimple 1978). The isotropic 
pressure tensor of the hot gas in early-type galaxies greatly 
simplifies measurement of the mass distribution over stellar 
dynamical methods. Moreover, since the shape of the vol- 
ume X-ray emission traces the shape of the gravitational 
potential independent of the (typically uncertain) gas tem- 
perature profile (Buote & Canizares 1994, 1996a), the shape 
of the mass distribution can be accurately measured in a way 
that is quite robust to the possible complicating effects of 
multi-phase cooling flows and magnetic fields (see Buote & 
Canizares 1997b). 

Presently, X-ray measurements of the mass distribu- 
tions in early-type galaxies are in hibited by limi tations in 
the available data. The ROSAT (|lrumpcr 1985|) Position 



Sen sitiv e Proportional Counter (PSPC) (Pfcffermann et al 



1987p h[a.s inadequate spatial resolution (PSF ~ 30" FWHM) 
to map the detailed mass distributions for all but the largest 
nearby galaxies, and the limited spectral resolution and 
band width complicates interpretation of the measured tem- 
perature profiles (Buote & Canizares 1994; Trinchieri et al. 
1994; Buote & Fabian 1997). Although equipped with su- 
perior spatial resolution (PSF ~ 4" FWHM), the ROSAT 
High Resolution Imager (HRI) (David et al. 1997) has too 



small an effective area and too large an internal background 
to provide images of sufficient quality for many galaxies for 
radii r J> R e . Among the few galaxies with detailed mea- 
surements of their radial mass profiles are NGC 507 (Kim & 
Fabbiano 1995), NGC 1399 (Rangarajan et al. 1995; Jones 
et al. 1997), NGC 4472 (Irwin & Sar azin 1996), NGC 4636 
(Trinchieri et al. 1994), NGC 4649 ( |Brighenti fc Mathews 
1997Qnd NGC 5044 (David et al. 1994). 

The shape of the gravitating mass has been measured 
via X-ray analysis for the E4 galaxy NGC 720 and the E7/S0 
galaxy NGC 1332 and found to be at least as elongated 
as the optical isophotes (Buote & Canizares 1994, 1996a, 
1997a). For NGC 720, which has more precise constraints, 
the ellipticity of the gravitating matter is e mass = 0.44 — 0.68 
(90% confidence) compared to the intensity weighted ellip- 
ticity of the optical light, (e) = 0.31 (Buote & Canizares 
1997a). In addition, the X-ray isophotes of NGC 720 twist 
from being aligned with the optical isophotes within R e to a 
position ~ 30° offset at larger radii. This twist, when com- 
bined with the ellipticities of the X-ray isophotes, cannot 



be explained by the projection of a reasonable triaxial mat- 
ter distribution and thus may implicate a dark matter halo 
misaligned from the stars (Buote & Canizares 1996b; Ro- 
manowsky & Kochanek 1997). 

NGC 720 and NGC 1332 were selected for analysis since 
they are isolated, significantly elongated in the optical, suf- 
ficiently bright, and sufficiently dominated by emission from 
hot gas in the ROSAT band. In this paper we present X-ray 
analysis of the classic "shell" galaxy, NGC 3923, which is 
the last galaxy of which we are aware that satisfies these se- 
lection criteria and has deep ROSAT observations. This iso- 
lated E4 galaxy has both archival ROSAT PSPC and HRI 
data and its ASCA spec trum has been analyzed previously 
( Buote fc Fabian 1997 ). This will serve as our final case 
study until the impending launch of AXAF revolutionizes 
this field. 

The organization of this paper is as follows. In ^ we 
describe the ROSAT observations and the data reduction. 
We discuss removal of point sources in ^ Measurements of 
the ellipticities of the X-ray isophotes and the radial pro- 
files are described in §^ and §^ respectively. Analysis of the 
PSPC spectrum is presented in §H. We give results for the 
Geometric Test for dark matter in §[j] and constraints on the 
shape and radial mass distribution from detailed hydrostatic 
models in Finally, in ^ we give our conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 PSPC 

NGC 3923 was observed with the PSPC for 14.5 ks from 
17-19 December 1991 and for 24 ks from 23-26 June 1993. 
Both observations were positioned at the field center where 
the point spread function is smallest. Since spatial resolution 
is of principal importance for our analysis of the shape of the 
X-ray surface brightness, we analyze only the PSPC data in 
the hard band (PI channels 42-201, E « 0.4 - 2.0 keV) to 
further optimize the size of the PSF; for details regarding 
the PSPC PSF s ee Ha singer et al. (1993, updated 1994) and 
see Aschenbach (1995) for a description of the ROSAT X-ray 
telescope. 

We reduced each observation separately with the stan- 
dard XSELECT, FTOOLS, and IRAF-PROS software. 
Firstly, the events files were cleaned of after-pulse signals by 
removing any events following within 0.35 ms of a precur- 
sor. We then removed large fluctuations in the light curves 
indicative of scattered light from the Bright Earth, Sun, or 
SAA; this resulted in filtered exposures of 13.5 ks for the 
1991 observation and 21.8 ks for the 1993 observation. 

To optimize signal-to-noise (S/N) and bin-size require- 
ments for computing ellipticities (see §k|), we binned the 
cleaned events files into images with 5 pixels. Exposure 
maps for the appropriate energy band (PI=42-201) were 
then generated for each image and were then used to create 
flattened images; note that the intrinsic resolution of the 
exposure maps is 15". Finally, we aligned the images us- 
ing bright point sources in the field and then combined the 
images. The final image is displayed in Figure |l| 
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2.2 HRI 

NGC 3923 was observed with the HRI for 45 ks from 10- 
17 January 1995 and for 26 ks from 22-26 June 1995. Both 
observations were positioned at the field center where the 
PSF is smallest (David et al. 1997). For detailed explanation 
of our reduction of HRI data we refer the reader to the 



related study of NGC 720 by Buote & Canizares (1996b) 



We restricted the data to those in pulse-height analyzer 
(PHA) bins 2-8 since they maximized the S/N of the data of 
each observation. For each observation we binned the events 
into images with 1" pixels to optimize determination of el- 
lipticities within a 30" radius of the galaxy center (§|i|). Ex- 
posure maps were generated with the standard software for 
each image and then used to flatten the images. 

An accurate aspect solution is critical to analysis of 
isophote shapes with the HRI on small scales (r 15"), with 
the amplitude of the asymmetrical distortion due to incom- 
plete aspect cor rection typically b eing most important for 
(r ~ 5" - 10") ( [David et al. 1997[ ). The S/N of each of the 
NGC 3923 observations are too small to usefu lly perform the 
aspect correction algorithm of Morse (1995). The low S/N 



also makes it difficult to make strong statements about the 
ellipticity of isolated point sources indicati ve of as pect error. 
As a result, following Buote & Canizares (1996b) we search 



for possible aspect errors by examining the positions of point 
sources for each of the OBIs; i.e. the time intervals during 
which the spacecraft continuously pointed on the target. 

Three OBIs contain most of the "on time" of the Jan- 
uary observation: 19 ks for OBI-6, 20 ks for OBI-11, and 9 
ks for OBI-14; 28 ks is distributed fairly evenly among 12 
other OBIs. Using OBI-11 as a reference, the bright point 
sources in the NGC 3923 field are displaced by 2" ± 1" for 
OBI-6, and 3" ± 1.5" for OBI-14. These shifts, though sta- 
tistically significant, are very consistent with the expected 
aspect uncertainties for a typical observation. To register all 
of the OBIs to one coordinate frame, we binned together 
those OBIs that were observed close together in time; i.e. 
OBIs 1-5, 6-9, 10-11, 12-13, and 14-15. Once registered to 



the same coordinate frame, the images were all added. (Note 
that the image for each of these OBI groups was binned as 
above.) 

The June observation proved to be problematic because 
none of the 19 OBIs had a long enough exposure to provide 
very accurate source positions. Hence, we were unable to 
perform a reliable test of the aspect solution. In order to 
make some estimate of the error, we grouped those OBIs 
closest together in time: OBIs 1-3, 4-12, 13-16, and 17-19. 
For three of these groups, less than 1" shifts were required, 
but the OBI 4-12 group required an 8" ± 1" shift. This large 
shift is questionable since all but one source was too faint 
to even obtain a centroid measurement using the alignment 
software in IRAF. 

As it is unclear from considerations of the OBIs how ac- 
curate is the aspect error for the June observation, we com- 
pared the ellipticity and orientations (computed as described 
in §^) of the surface brightness with those of the January 
observation. Unfortunately, we find significant disagreement 
between the two observations; note the disagreement occurs 
whether or not we include the OBIs 4-12 for the June obser- 
vation. In particular, the position angles determined from 
the January observation are fully consistent with the PSPC 
data and the optical position angles for r ^ 30", which 
agrees with NGC 720 and NGC 1332 (Buote & Canizares 
1996a,b); the ellipticities are also consistent with the PSPC 
data. The June observations, however, have position angles 
that differ by over 60° and the ellipticities, particularly for 
r J; 20" — 30" are less than 0.1 as opposed to ~ 0.25 for the 
January observation. Furthermore, the radial profile (see §^|) 
of the June observation is significantly flatter for r £ 10". 

Thus, the isophote shapes and orientations are incon- 
sistent for the two observations. The differences point to a 
serious aspect error in the June observations. That is, al- 
though we would expect aspect error to induce ellipticity in 
an intrinsically circular source, if the distortion occurs nearly 
along the minor axis of a moderately elliptical source, it will 
reduce the ellipticity of the source, as well as alter the posi- 
tion angle. The flatter inner radial profile of the June data is 
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Table 1. 


Identified Point Sources 


Source 


R. A. 


Decl. 


1 


ll h 51 m 06 s 


-28° 46' 47 


2 


11 50 58 


-28 43 58 


3 


11 50 15 


-28 53 22 


4 


11 51 04 


-28 57 30 


5 


11 51 29 


-28 43 51 


6 


11 51 33 


-28 47 41 


7 


11 51 35 


-28 56 23 


8 


11 51 40 


-28 46 24 


NWa 


11 50 53 


-28 47 04 


NWb 


11 50 58 


-28 46 30 



These sources (expressed in J2000 coordinates) were identified 
from visual inspection of the PSPC image within a 10' radius of 
the galaxy center. They were removed from the image for spatial 
analysis as described in ^ 

also consistent with uncorrected wobble. As a result, we do 
not include the June data in our analysis. (This is not critical 
since the S/N to be gained over the 45 ks is marginal.) 

We display the reduced January observation in Figure 



sumed symmetry of the hot gas distribution. If the gas is 
approximately ellipsoidal, then we can replace source (1) 
with the corresponding emission obtained by reflecting the 
source over the galactic center; i.e. essentially on the other 
side of the major axis (a — > —a). Fortunately, we find that 
the ellipticities and position angles are virtually the same 
whether we remove source (1) by subtracting a model for 
the local background or by symmetric substitution. (We use 
the former method for ensuing analysis.) 

The emission associated with NWa,b cannot be so eas- 
ily removed because it is extended, and it is not obvious 
how to define a background model. Fortunately this emis- 
sion only begins to affect the ellipticities for semi-major axes 
<; 120". At these distance the S/N does not allow accu- 
rate constraints. However, this emission does need to be re- 
moved from the radial profile, and thus we iterate the local 
background method to remove the emission associated with 
NWa,b. This method is suitable for the radial profile since 
the azimuthal averaging is not overly sensitive to small non- 
axisymmetric residuals. (We mention the effect of removing 
this emission on the radial profile in §@.) 



3 POINT SOURCES 



It is readily apparent from inspection of Figure particu- 
larly for the higher S/N PSPC data, that the X-ray emis- 
sion of NGC 3923 is significantly contaminated by fore- 
ground/background point sources. Within a 10' radius of 
the galaxy center at least 8 point sources are easily identi- 
fied by visual inspection of the PSPC image. (Only r ^ 3' is 
shown in the Figure.) We list the positions of these sources 
in Table [l[ 

For our analysis of the isophote ellipticities and orien- 
tations (§^) and the radial profile (§^) of the X-ray surface 
brightness, we wish to analyze only the distribution of the 
diffuse emission associated with NGC 3923 and thus these 
contaminating sources must be removed. The most impor- 
tant contaminating source is number (1) (which lies along 
the optical major axis to the N-E) since it lies closest to the 
galactic center where the S/N , and thus constraints on the 
surface brightness, are best; this source is also apparent in 
the HRI image. Next in importance is the extended emis- 
sion approximately 2.5' to the N-W along the optical minor 
axis which appears to consist of emission from at least two 
point sources. These candidate sources are listed as NWa,b 
in Table |l[ 

Our preferred method for removing sources, which is 
well s uited for analyzing qu adrupole moments of X-ray im- 
ages (Buote & Tsai 1996), is to first choose an annulus 



around each source to estimate the local background. Then a 
second order polynomial surface is fitted to the background 
which then replaces the source. We removed sources 2-8 us- 
ing this method, but sources (1) and NWa,b require some 
elaboration. 

Since source (1) affects ellipticity measurements for 
r <; 100" we examined how robust the ellipticity and po- 
sition angle were to the method used to remove the source. 
Another method to remove sources is by "symmetric substi- 
tution" (see Strimple & Binney 1979; Buote & Tsai 1995; 
Buote & Canizares 1996a). This method exploits the as- 



4 X-RAY ISOPHOTE SHAPES AND 
ORIENTATIONS 

As is typical for current X-ray data of early-type galaxies 
the small number of counts (~ 1000) for the PSPC and 
HRI images of NGC 3923 implies that we can only hope 
to measure with any precision the ellipticity and position 
angle of the aggregate X-ray surface brightness in a large 
aperture. The method we employ is an iterative procedure 
(| Carter fc Metcalfe 1980| ) and is analogous to computing 



the two-dimensional moments of inertia within an elliptical 
region where the ellipticity, cm, is given by the square root of 
the ratio of the principal moments and the orientation of the 
principal moments gives the position angle, 9 m (see Buote 
& Canizares 1994 for application to ROSAT images). The 
parameters €m and 9 m are good estimates of the ellipticity 
(e) and position angle (9, or P.A.) of an intrinsic elliptical 
distribution of constant shape and orientation. For a more 
complex distribution 6m and 9m are average values weighted 
heavily by the outer parts of the regions. 

We estimate the uncertainties on 6m and 9m using 
a Mont e Carlo procedure described in Buote & Canizares 
(1996b). In sum, the procedure involves constructing 1000 



realizations of the PSPC and HRI images taking into ac- 
count statistical noise and unresolved sources; these unre- 
solved sources are modeled according to th e log N(> S) — 
log 5* distribution given by Hasinger et al. (1993) and their 
profiles are given by the appropriate PSPC or HRI PSF. 
The 90% confidence limits on 6m, for example, are denned 
by the 5th and 95th percentile values computed from the 
1000 simulations. 

The profiles of em and 9m are listed for the PSPC in 
Table § and for the HRI in Table |. We also display the e M 
profiles and their 68% uncertainties in Figure 0. The aper- 
tures listed in the tables are chosen so that each increment 
in semi-major axis a consists of approximately 100 source 
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Table 2. PSPC Ellipticities and Position Angles 
o 

(arcsec) e M 68% 90% 6 M 68% 90% Counts 



30 0.10 0.04-0.19 0.03-0.22 42 19-55 11-77 1116 

40 0.12 0.09-0.15 0.07-0.17 25 15-41 07-53 1331 

50 0.07 0.05-0.12 0.03-0.13 31 20-52 09-66 1460 

60 0.08 0.05-0.15 0.03-0.19 57 38-81 26-101 1564 

75 0.15 0.12-0.19 0.09-0.21 51 42-58 36-65 1670 

90 0.15 0.11-0.20 0.06-0.23 48 37-60 31-71 1754 

110 0.10 0.08-0.16 0.05-0.19 39 26-61 13-156 1870 



The values of tyi (and confidence limits) are computed within an aperture of semi- 
major axis a on the image with the background included; the counts, however, 
have the background subtracted. The values of 9m are given in degrees N through 
E. 



Table 3. HRI Ellipticities and Position Angles 
a 

(arcsec) e M 68% 90% M 68% 90% Counts 



12 0.25 0.13-0.32 0.08-0.39 69 44-83 30-92 309 

17 0.26 0.14-0.31 0.08-0.36 59 50-76 43-88 412 

23 0.22 0.09-0.27 0.05-0.32 66 47-81 34-102 506 

32 0.13 0.07-0.21 0.04-0.26 63 35-82 15-107 631 

44 0.14 0.11-0.27 0.07-0.31 76 41-87 31-103 735 

60 0.16 0.10-0.28 0.05-0.33 74 44-101 33-123 825 



See Table [| 



counts. We restrict the PSPC profile to a < 110" because of 
contamination from sources NWa,b (see §^|)Q. 

The values of em for the HRI have considerably larger 
estimated uncertainties than for the PSPC. In their regions 
of overlap, the HRI and PSPC give consistent ellipticities. 
For the smallest a ^ 20", the HRI gives some indication 
of large ellipticity, £m <; 0.20, but the 90% lower limits are 
less than 0.10; at these radii residual aspect errors could 
be important. At these a the PSPC is blurred by the PSF, 
but at a ~ 70" - 90" (~ 10 kpc) the PSPC gives the best 
constrained ellipticity for either the HRI or PSPC of em ~ 
0.15 ± 0.05. Although significant, this em is less than that 
measured on similar scales for NGC 720 (~ 0.25; Buote & 
Canizares 1994; 1996b) and for NGC 1332 (~ 0.20; Buote 
& Canizares 1996a). 

Like NGC 1332, the P.A. profiles of the PSPC a nd HRI 
data ar e con sistent with the optical value of ~ 49° (Jedrze- 
jewski 1987), albeit within the rather large uncertainties. 
Moreover, the PSPC P.A. for a <; 90" may be affected by 
the bright source (1) and NWa,b (see Table g), and thus it 
would be premature to rule out a P.A. twist like that in NGC 
720 until higher resolution, better S/N data (e.g. AXAF) is 
obtained. 

Finally, over the region a < 110" we find no obvious 



t We emphasize that the bright point source (1) (see Table |lj 
and Section ^) does not affect the ellipticity measurements for 
a < 90"; i.e. we also measure em = 0.15 for a = 75", 90" on the 
image including source (1). However, the values at larger radii 
differ considerably depending on whether the source is removed 
or not; e.g. at a = 120", we obtain em = 0.3 with the source and 
e M = 0.1 without the source. As a result, when comparing model 
ellipticities to data in future sections of this paper we use only 
the ellipticities for a < 90" which are robust to the inclusion or 
exclusion of source (1). 



asymmetries in the (PSPC or HRI) surface brightness (Fig- 
ure ji]) suggestive of unresolved sources or environmental ef- 
fects like ram pressure or tidal distortions; e.g. the centroids 
of the apertures are very steady over the semi-major axis 
listed in Tables and g. We also have divided up the re- 
gion r < 75" into quadrants (aligned with the major axis of 
NGC 3923) and have found the counts in each quadrant to 
be consistent with each other within (1 — 2)cr errors, similar 
to the agreement seen for NGC 720 (Table 2 of Buote & 
Canizares 1996b). 



5 RADIAL PROFILE OF X-RAY SURFACE 
BRIGHTNESS 

In order to construct the azimuthally averaged radial pro- 
files for the PSPC and HRI data we require measurements 
of their respective background levels. We selected annular 
regions centered on the galaxy that are sufficiently far from 
the galaxy center so that contamination from the galaxy is 
minimal. (Also, any point sources were removed.) We obtain 
a background rate of 2.4 x 10 -4 cnts s _1 arcmin -2 for the 
PSPC and a rate of 3.3 x 10" 3 cnts s" 1 arcmin for the 
HRI. 

We binned the radial profiles so that each bin had ap- 
proximately the same S/N. However, for the innermost bins 
the S/N is larger because we did not oversample the respec- 
tive PSFs. The centers of the radial bins for both the PSPC 
and HRI data were determined by the centroid of the circu- 
lar region containing ~ 80% of the total flux; in neither case 
did this choice critically affect the radial profile shape. We 
display the radial profiles in Figure [| along the the PSFs; 
note that residuals of the emission near sources NWa,b may 
be seen in the PSPC profile for r ~ 130" - 150". The X-ray 
emission is clearly extended for both data sets. 



© 0000 RAS, MNRAS 000, 000-000 



6 D. A. Buote and C. R. Canizares 



-i — i — i — i — r- 



i 1 1 1 i 



R-band 



PSPC 



4> 



_i_ 



20 40 60 80 

a (arcsec) 



100 



120 




40 60 80 

a (arcsec) 



120 



Figure 2. Moment cllipticities and lcr errors for PSPC (left) and HRI (right) data. The isophotal ellipticities of the .R-band light 
from Jedrzejewski (1987) are indicated by the dotted line. The dashed line represents the isopotential ellipticities of, <E>ft, the potential 
assuming M oc Lr (see §[iT). 
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Figure 3. Radial surface brightness profiles of the PSPC (left) and HRI (right) data. The error bars indicate the data and the horizontal 
lines through each point show the bin sizes. The solid horizontal lines across each plot indicate the background level. The best-fit /3 
model (filled circles), the PSFs (crosses), and the .R-band light convolved with the respective PSFs (boxes) are also shown. All of these 
quantities have been binned as the X-ray data. 



A convenient parametrisation of the X-ray radial pro- 
files of early-type g alaxies is given by the "/3 model" , oc 
{Rl + i? 2 )-3/3+o.5 ( [Cavaliere fc Fusco-Femiano 1976| ). The 



model assumes that the gas is isothermal and that the stars 
considered as test particles, follow a King law. Although 
these assumptions are unlikely to be strictly valid for el- 
lipticals, this model provides a reasonabl e description of the 
ROSAT radial profiles of both NGC 720 (|Buote fc Canizares 
1996b[]and NGC 1332 ( [Buote fc Canizares 1996a| ) as well 
as many other galaxies (e.g. Forman et al. 1985). Moreover, 
for ROSAT data the /3 model typically gives fits of quality 
very similar to more sophisticated mass models and thus it 
serves as a good benchmark for these more general models 
(as in §§). 



Table 4. Simple j3 Model Fits to Radial Profile 



Re 



Data 


(arcsec) 


P 


x 2 


dof 


^red 


PSPC 


"5 fi +11 
°-°-l.O 


n 47+0.02 
u - £ " -0.02 


8.9 


9 


1.0 


HRI 


tr 1+1.6 
0.1_ 1 2 


46+ - 05 


1.4 


1 


0.4 


PSPC + HRI 


5 4+ 08 


46+ - 02 
u -^ D -o.oi 


10.7 


15 


0.7 



The best-fit values and 90% confidence limits on one interesting 
parameter are listed for R c and j3. The PSFs of each data set 
have been incorporated into the fits. 

In Table [| we list the results of fitting the fi model 
(convolved with the appropriate PSF) to the X-ray data; 
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Figure 4. Absorbed single-temperature MEKAL model fit jointly to the inner (upper curve) and outer (lower curve) regions (see 
Only the June, 1993 data is shown. 



only those bins with S/N > 5 were included in the fits. 
When fitted separately to the PSPC and HRI data, the j3 
model gives a good quality fit and the derived R c and j3 
parameters are similar; note that if the emission from sources 
NWa,b is not removed from the PSPC data, the parameters 
of these fits are essentially unaffected but the value of x 2 
increases to 17.4. 

The best-fit model obtained from jointly fitting the 
PSPC and HRI data is shown in Figure |[ The {3 model 
provides a reasonably good fit over a large radius range for 
both data sets. For com parison, we plot the ii-band data 
from Jedrzejewski (1987) in Figure H. We have binned the 



optical data like the X-rays and have convolved with the ap- 
propriate PSFs; for radii larger than the limiting a — 103" of 
Jedrzejewski, we extrapolate the data using the best-fitting 
De Vaucouleurs model. For radii larger than the PSFs, the 
shape of the optical profile, E_r, is very similar to the shape 
of Ti x in agreement with expectations from steady-state cool- 
ing flow models (e.g. Sarazin 1987). However, the X-ray and 
optical profiles disagree markedly for small radii (r ^ 10") 
with E_r being much more centrally peaked than T, x . Be- 
cause of its smaller PSF this effect is more pronounced in 
the HRI data. 

In particular, it is clear that a population of discrete 
X-ray sources distributed like E_r cannot contribute signifi- 
cantly to the emission. If we assume that the hot gas is de- 
scribed by a ft model or related model (see §[]) then adding 
such a discrete component only worsens joint fits to the 
PSPC and HRI data. We find that f h Jf disc > 3.6 (90% 
confidence), where fhg is the flux of the hot-gas component 
(/3 model) and fdisc is the flux of the discrete sources where 
the 0.4-2 keV flux is computed within a circle of r = 2'. 

Spect ral an alysis of ASCA data of NGC 3923 by Buote 
& Fabian (1997) shows that two temperature models are re- 
quired by the data with a cold component, Tc = 0.55 keV, 
and a hot component, T H = 4.2(> 2.2) keV (90% confi- 



dence); the ratio of the 0.5-2 keV flux of the cold component 
to the hot component is 1.9 (1.3-2.8) at 90% confidence. If 
the discrete sources are indeed distributed like E_r, then the 
emission of the hot component cannot be entirely due to dis- 
crete sources. This would indicate that the emission of the 
hot component is actually largely due to another phase of 
the hot gas and that Th is overestimated due to an artifact of 
fitting low S/N data as suggested by Buote & Fabian. (Note 
that the cold-to-hot flux ratio derived by Buote & Fabian 
in the 0.5-2 keV band remains essentially the same when 
computed in the 0.4-2 keV band analyzed in our paper.) 



6 SPECTRAL ANALYSIS 

Spectral analysis of the X-ray data is required for our study 
of the mass distribution in NGC 3923 to determine (1) how 
much of the emission is due to hot gas, and (2) the temper- 
ature profile of the gas. The superior spectral resolution of 
ASCA is better suited than ROS AT to address issue (1), and 
as discussed at the end of the previous section, the ASCA 
spectrum is consistent with ~ 35% of the X-ray emission 
in the 0. 5-2 keV band arising fr om a population of discrete 
sources (Buote & Fabian 1997). Since this much discrete 



emission is inconsistent with the ROS AT radial profiles we 
shall neglect it and assume all of the emission arises from 
hot gas. This assumption does not seriously affect analysis 
of the mass distribution if the discrete contribution is ^ 20% 



(Buote & Canizares 1997a) 



Unlike the ASCA data, the ROSAT PSPC data allows 
us to directly measure any temperature gradients. To in- 
vestigate this issue we analyzed the PSPC spectra in a 
circular region with r — 30" and an annular region with 
r — 60" — 120"; the spectra of these regions are shown in 
Figure ^. We fit a model consisting of Galactic absorption 
(N H = 6.4 x 10 20 cm" 2 - Stark et al. 1992) and a thin 
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Table 5. Spectral Fits 



Name 


T 


z 


x 2 


dof 


y2 




(keV) 










0" - 30" 


n cn+0.05 
u - ou _0.04 


0- 44 to.?5 


51.3 


50 


1.0 


60" - 120" 


^S+ 017 
U - M _0.15 


09+ 67 


7.9 


12 


0.6 


BOTH 


50+ - 04 
U - OU -0.05 


45+1-2 


61.6 


61 


1.0 



Results of fitting a MEKAL model modified by Galactic photo- 
electric absorption. We list 90% confidence limits on one interest- 
ing parameter. 

t hernia l plasma. Our thermal plasma model is the MEKAL 
model ]vhich is a modification of the original MEKA code 



(Mewe, Gronenschild, & van den Oord 1985; Kaastra & 
Mewe 1993) where the Fe-L shell transitions crucia l to the 
X-ray e mis sion of ellipticals have been re-calculated (Liedahl 
et. al 1995). We ta ke solar abundances according to Anders 
& Grevesse (1989) and photo-electric absorption cro ss sec - 
tions according to Balucihska-Church & McCammon (1992). 



All of the spectral fi tting was imp lemented with the 
software package XSPEC (Arnaud 199q) using the \ 2 mini- 



mization method. In order for the weights to be valid for the 
X 2 method we regrouped the PI bins such that each group 
contained at least 20 counts. We restricted the fits to en- 
ergies above 0.2 keV because of the low S/N coupled with 
uncertainties in the PSPC response at these lowest energies. 
We fitted the December and June observations jointly, but 
with their normalizations as free parameters. 

The results of the spectral fits are given in Table ^. The 
model is an excellent fit to the low resolution PSPC data 
in both regions with the temperatures and abundances giv- 
ing consistent values within the uncertainties. There is some 
indication that the inner region prefers higher abundances 
than the outer region. This may suggest a real abundance 



gradien t, or it may be due to a difference in S/N ( Buote & 
Fabian 



1997) 



In any event, a joint fit to the inner and outer 
region is formally acceptable and gives parameters similar 
to the inner region. Hence, we find no evidence of tempera- 
ture gradients, and since the PSPC spectral constraints are 
even tighter than found for NGC 720 any gradients consis- 
tent with this data are unimportant for analysis of the mass 
distributions (see Buote & Canizares 1994). 

Note that we did jointly fit the PSPC and ASCA data 
and found that the constraints did not change apprecia- 
bly from the ASCA results alone, even for two-temperature 
models. 



7 GEOMETRIC TEST FOR DARK MATTER 

The shapes of the X-ray isophotes allow the shape and radial 
distribution of gravitating matter to be probed in a way that 
is more robust than the traditional spherical approach (for 
a review see Buote & Canizares 1997b). Assuming the hot 
gas is approximately in hydrostatic equilibrium and that the 
emission is adequately described by a single phase, then the 
volume X-ray emissivity, j x , obeys an "X-ray Shape The- 
orem" (Buote & Canizares 1994, §3.1; 1996a, §5.1) which 
states that j x and the gravitational potential, have the 
same three-dimensional shapes independent of the temper- 
ature profile of the gas. One may thus make a "Geometric 



Test" for dark matter that is distributed differently from 
L (with L the luminosity distribution of the optical stellar 
light) by comparing the shape of $_l generated assuming 
M oc L with the shape of j x obtained from deprojecting the 
X-ray image data. As this comparison is independent of the 
poorly constrained temperature profile, T(r), of the gas (see 
§^) this test for dark matter is more robust than computing 
the radial mass distribution which is directly proportional 
to T(r)|. 

NGC 3923 has negligible stellar rotation, but if rotation 
of the gas is dynamically important then we must replace <£> 
with the appropriate effective potential. Theoretically, even 
witho ut strong st ellar rotation one may (Kley & Mathews 



1995) or may not (Nulsen et al. 1984) expect a rotating cool- 



ing flow to develop depending on whether angular momen- 
tum of the gas is conserved. Highly flattened X-ray isophotes 
indicative of a cooling disk have not been observed in ellip- 
ticals and thus we shall ignore any contribution from gas 
rotation^]. 

To apply the Geometric Test we constructed a constant 
M/L potentia l usin g the J?-band surface photometry, E_h, of 
Jedrzejewski (1987). The major-axis profile of En is well fit- 
ted by a De Vaucouleurs Law with an effective semi-major 
axis, a e — 92". Taking into account the ellipticity of the 
isophotes, we fitted the De Vaucouleurs Law to the mean 
radial profile, where the mean radius is, r = Vab = a^fq, 
where a is the major axis and q the axial ratio. This gives a 
mean effective radius, R e = 73". Since the 3-D density giv- 
ing rise to a De Vaucouleurs Law is well approximated b y 
the Hernquist profile, p oc r~ 1 (r s + r)~ 3 ( Hernquist 1990 ), 
we consider for our M oc Lr model a Hernquist profile 
with r s = i? e /1.8153 = 40". Furthermore, we take p to 
be stratified on oblate or prolate spheroids of constant ellip- 
ticity, where we use the intensity weighted ellipticity of E_r, 
(e R ) = 0.30. 

This simple prescription for the M oc Lr model is suf- 
ficient for comparison to the ROSAT X-ray data because of 
the relatively crude X-ray constraints: i.e. only a global con- 
straint on the isophote shapes of the PSPC data over a small 
range in radius, r ~ 70" — 90", is useful for the comparison 
- as we show below, the HRI data do not provide inter- 
esting constraints over the PSPC data. With better quality 
data from future missions which accurately measure ellip- 



t The X-ray shape analysis is also more robust to the effects 
of mult i-phase gas and magnetic fields - see Buote & Canizares 
rtl997b|). 



° We should also emphasize that the wiggles displayed by the 
smoothed X-ray isophotes in Figure [l] are not statistically sig- 
nificant and thus do not suggest departures of the gas from hy- 
drostatic equilibrium; i.e. the S/N of the PSPC and HRI data 
is too poor to measure deviations of isophote shapes from ellip- 
tical. However, even if those wiggles in the contours were due 
to real non-equilibrium motions in the gas, the hydrostat ic as - 
sumption is still very likely a good one. Buote &; Tsai ( 1995[ ) 
analyzed the X-ray isophotes of an X-ray cluster formed in an N- 
body/hydrodynamical simulation. They concluded that relatively 
soon after major mergers the cluster settled down quickly to a re- 
laxed state which allowed accurate determination of the shape of 
the mass from analysis of the X-ray isophote shapes, even though 
the simulated clusters had X-ray isophotes far from elliptical in 
shape. 
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Figure 5. Geometric Test for Dark Matter 
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ticity gradients over a large range in radius a more accurate 
M oc Lr approximation will need to be considered. 

We show in Figure || the ellipticities of the isopoten- 
tial surfaces of the potential, <&r, assuming M oc Lr and 
edge-on oblate symmetry. These ellipticities are considerably 
smaller than those of the _R-band light because the spheri- 
cally symmetric monopole term in the potential dominates 
the isopotential ellipticities for r <; r c p[ The ellipticities 
of the X-ray surface brightness, particularly for the PSPC 
data for a ~ 70" — 90", appear to be significantly larger 
than the isopotential ellipticities and thus indicate a failure 
of the M oc Lr hypothesis. 

However, to rigorously compare <&r to the X-ray data 
we must formally deproject E^ to get the ellipticities of j x . 
The procedure we adopt for this comparison, which is suit- 
able for the relatively crude constraints provided by the X- 
ray data, is to first fit a simple model to the radial profile of 
E x ; for this purpose, the results of the /3 model fit jointly to 
the PSPC and HRI data are suitable (see §H). This model 

for T, x (R) is then deprojected to obtain j x (r|fl. We then as- 
sign to j x the ellipticities of $u which gives a spheroidal 
emissivity distribution, j x (r, 9) . This spheroidal j x is then 
projected back onto the sky plane, while adjusting R c and (3 
to maintain a best fit of the radial profile of S x (convolved 
with the appropriate instrument PSF), to yield a M oc Lr 
model of the X-ray surface brightness. We compute moment 
ellipticities, € x (&r), for 1000 Monte Carlo simulations of this 
model in analogy to the data (see §Q). 

The la error bars for e x (&R) as well as the observed X- 
ray ellipticities and the 3-D isopotential ellipticities of "1>_r 
are displayed in Figure ^. Although most of the measured X- 
ray ellipticities lie within the error bars, the best measured 
PSPC ellipticities from a ~ 70" - 90" exceed the values 
predicted by the edge-on oblate $r at the 85% confidence 



level. The edge-on prolate $_r is formally discrepant at the 
80% level at these radiQ 

Hence, the X-ray isophote shapes indicate that the 
e x (&R) model is either too round, too centrally concen- 
trated, or both; i.e. dark matter is required which is flat- 
tened and probably more extended than E.r. This discrep- 
ancy, though formally marginal and of lesser significance 
than found for the other two early-type galaxies studied 
NGC 720 and NGC 1332 (e.g. Buote & Canizares 1997a,b), 
is of precisely the same character: i.e. flattened and extended 
dark matter is required in these ellipticals. Since the flattest 
optical isophotes of NGC 3923 are rounder than those of 
NGC 720 and NGC 1332, it is possible that the symmetry 
axis of NGC 3923 is inclined more along the line-of-sight 
than the other two. This would account for the marginal 



** We caution that interpretation of the significance of this dis- 
crepancy must take into account the PSPC PSF as we have done; 
i.e. the two deviant data points for a ~ 75" — 100" are not sim- 
ply random fluctuations weighted equally with other points over 
the radial range investigated. (Also, care is required in the in- 
terpretation of the last measured ejvf at a = 110" for the PSPC 
data since its value depends on the bright point source that is 
removed — see Sections jj] and ^.) This point is illus trated in Fig- 
ures 13 (a) and (b) of Buote & Canizares (1996a) for a similar 



1 See Figure 2-13 in §2.3 of Binney & Tremaine ( 1987 ) for a 
related discussion of isopotential ellipticities. 
I Note that the parameters in Table Q correspond to the depro- 
jected parameters. 



X-ray analysis of NGC 1332. That is, the X-ray ellipticity profiles 
produced by different models are smeared out by the PSPC PSF, 
and thus a discrepancy is only achieved at radii large enough so 
that the ellipticity differences between models exceeds the rela- 
tively large error bars of the measured ellipticity; of course, the 
radial range is limited from above by the decreasing S/N. Thus, 
the radii a ~ 75" — 100" are where the different models con- 
volved with the PSPC PSF begin to show significant ellipticity 
differences and the X-ray data still have good constraints on the 
measured ellipticity. NGC 3923, NGC 1332, and NGC 720 are 
quite similar in this regard since they have similar length scales 
and similar S/N PSPC observations. However, the HRI data of 
NGC 720 are of sufficient quality so that the effect of the PSFs 
on the Geom etric Te st can be clearly seen: see Figure 2 of Buote 
& Canizares (1997b). 
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significance of the NGC 3923 result if the intrinsic shapes of 
the three galaxies are in fact similar. 

We no te that , as explained in Section 5.1 of Buote & 
Canizares (1996a), if an ellipsoidal galaxy is inclined along 
the line of sight, then the Geometric Test results give a lower 
limit to the true discrepancy. That is, if mass follows light, 
then j x and <E>l must be co-axial and thus share the same 
inclination angle, i. Since the projection of an ellipsoid with 
i < 90° is necessarily rounder than if the galaxy were viewed 
edge-on, deprojection of E x and S_l assuming i = 90° will 
yield j x and which are rounder than in reality. However, 
this edge-on deprojection of an inclined ellipsoidal galaxy 
only means that any differences in the inferred shapes of j x 
and are smaller than if the true inclination angle were 
used. 



8 DETAILED HYDROSTATIC MODELS 

The results of the Geometric Test from the previous sec- 
tion suggest the presence of a dark matter halo that is flat- 
tened and more extended than the optical light. To generally 
constrain the allowed distribution of gravitating matter, we 
must employ explicit solutions of the hydrostatic equation. 
As the low S/N of the PSPC and HRI images does not jus- 
tify a sophisticated non-parametric inversion of the data, 
we instead consider simple intuitive models to place con- 
straints on the aggregate shape and radial distribution of 
gravitating mass. As in our previous studies (e.g. Buote & 
Canizares 1997b), we f ollow the pioneering approach of Bin- 
ney & Strimple (1978) and solve the hydrostatic equation 
assuming a single-phase, non-rotating, ideal gas, 



Pa 



exp [r (l - 4) 



(1) 



where T = um p &(0) /ksT and where, e.g. p g = p g (x)/p g (0). 
Equatio n (M[) assumes the gas is isothermal, but Strimple & 
Binney (1979) and others (Fabricant, Rybicki, & Gorenstein 
1984; Buote & Canizares 1994; Buote & Tsai 1995) have 
shown that the constraints on the shape of the gravitating 
matter distribution are not very sensitive to temperature 
gradients. Hence, we shall assume an isothermal gas which 
is consistent with the spectral constraints for NGC 3923 (§^|) 
and the temperature profiles of other galaxies (e.g. Buote & 
Canizares 1994; for a review see Sarazin 1997). 

Our procedure begins with a model for the gravitating 
mass. We consider spheroidal density distributions whose 
isodensity surfaces are concentric, similar spheroids. By ex- 
amining oblate and prolate configurations we bracket the 
intermediate behavior of triaxial models; i.e. as our impor- 
tant constraints on the models are the X-ray isophote shapes 
for semi-major axes a ~ 75", the detailed radial behavior of 
the isophote shapes and orientations which distinguish tri- 
axial models cannot be usefully constrained by the NGC 
3923 PSPC and HRI data. 

Once the type of density model and its associated ellip- 
ticity, e ma3S , are chosen, we compute $ and p g . (We initially 
focus on one-component models of the density of the grav- 
itating matter and then consider separate density models 
for the dark and luminous matter.) The model X-ray sur- 
face brightness is then generated by integrating p 2 g along the 



line of sightp| assuming the galaxy is viewed edge-on^]. We 
model the radial mass as either a softened isothermal po- 
tential with mass density, p oc (a 2 + a 2 ) -1 , or a Hernquist 
profile, p oc a~ 1 (a s + a)~ 3 , where a s is the appropriate mass 
scale length in each case. These models bracket the interest- 
ing range of densities in our previous X-ray analyses of NGC 
720 and NGC 13320. We convolve the surface brightness 
with the PSF of the appropriate detector and evaluate the 
model on a grid of pixels identical to that used for analyz- 
ing the data (§^). We construct the radial profiles as for the 
data (§^) and determine a 3 and T by fitting the model pro- 
file jointly to the PSPC and HRI data. (The normalizations 
of the PSPC and HRI are free parameters.) By comparing 
the ellipticity of the model surface brightness to the data we 
constrain the input ellipticity of the gravitating mass. 

In Table |^ we give the constraints on the shape of the 
gravitating matter obtained from these models. The mass is 
significantly flattened, e masa w 0.5, with the derived shapes 
being essentially the same for both the p ~ r~ 2 and Hern- 
quist models. These ellipticities are inconsistent with the 
intensity- weighted ellipticity of the i?-band light ((ej?) = 
0.30) at more than the 90% level and are marginally incon- 
sistent with the flattest optical isophotes ((e^ ax ) = 0.39). 
The derived range of values of e, nass and the fact that they 
exceed the average ellipticity of the optical isophotes are 
consistent with the previous results we obtained for NGC 
720 (and NGC 1332) (e.g. Buote & Canizares 1997a,b). 

In Figure ^ we plot the radial profiles of typical p ~ r~ 2 
and Hernquist models. The p ~ r~ 2 model fits are very sim- 
ilar in quality to those of the (3 model (§^) but with slightly 
larger values of x 2 - (The differences are not obvious from 
visual inspection.) However, the Hernquist fits are notice- 
ably worse than the p ~ r~ 2 model and have unacceptable 
Xred ~ 1-9- For r 15" the Hernquist model is too flat and 
then too steep out to r ~ 100". At larger radii, the small 
fitted values of F ~ 5.8 force the Hernquist radial profile to 
be flatter than the data. 

These deviations of the Hernquist model are most pro- 
nounced for the M oc Lr model (see §Q) which has a a = 40" 
and a best-fit F = 5.5 and \ 2 = 56.0 (16 dof). The PSPC 
ellipticities of the M oc L R model are 0.07 for a ~ 70" - 90" 
which is below the 90% confidence limits of the data (Ta- 
ble ^). Hence, by assuming the gas is isothermal we find 
that the M oc Lr model is inconsistent with the data at a 
significance level greater than indicated by the more robust 
Geometric Test (§|). 

It should be emphasized that most of the disagreement 
between the p ~ r~ 2 and Hernquist models occurs in the in- 

tt The plasma emissivity is a very weak function of temperature 
when convolved with the spectral response of either the PSPC or 
HRI. 

t+ We assume the symmetry axis of the galaxy spheroid lies in 
the plane of the sky; i.e. we are not attempting in this analy- 
sis to uncov er the true three-dimensional shape, though Binney 
& Strimple (|l97^) have found that the X-ray analysis is not ex- 
tremely sensitive to small inclination of the sym metry axis with 
respect to the sky plane. See Buote & Tsai (1995) for a thorough 
discussion of projection effects on X-ray shape analysis. 

The Hernquist density gives fits to the density profiles of halos 
in Cold Dark Matter simulations that are very similar to the 
universal model of Navarro et al. (1997). 
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Table 6. Ellipticity of the Gravitating Matter 



Model 
p ~ r 
Hernquist 



Oblate e mass 
68% 90% 
0.43-0.59 0.35-0.65 
0.45-0.61 0.36-0.66 



Prolate <= mass a s 

68% 90% (arcsec) |r| \ 2 dof 

0.40-0.54 0.33-0.58 2.2-3.1 6.6-7.2 16 15 

0.42-0.55 0.33-0.59 53-75 5.7-5.8 28 15 



Y 2 
1.1 

1.9 



Derived shapes of the gravitating matter for the spheroidal mass models assuming an isothermal gas (see 
§^). 90% confidence ranges are quoted for the scale length a s and |T|. Typical values of \ 2 are listed in 
the 90% intervals. Note for the M cc Lr model we obtain \ 2 = 56.0 for 16 dof. 
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Figure 6. Radial surface brightness profiles of selected oblate models fit to the PSPC (left) and HRI (right) data. The error bars 
indicate the data and the horizontal lines show the bin sizes. The filled circles represent the best-fit p ~ r -2 model having e ma ss = 0.50; 
the boxes represent the best-fit Hernquist model with e maS s = 0.50; the dotted line is the best-fitting constant M/L model. The radial 
profiles of these models have been binned as the X-ray data. 



nermost bins. For example, if we exclude the two inner bins 
for the PSPC (r = 0" - 30") and the three inner bins for the 

HRI (r = O"-15")0 we obtain a minimum X 2 of ~ 10.0 for 
the p ~ r~ 2 model and \ 2 ~ 11A for the Hernquist model. 
Since in these innermost regions the surface brightness pro- 
file could be modified by a multi-phase cooling flow or by 
magnetic fields, we suggest some caution in interpreting \ 2 
values of our models in these central regions until better 
data can distinguish between different scenarios. 

However, there is reason for optimism on both of these 
accounts. The distortion due to a multi-phase cooling flow 
is likely to be insignificant since single-phase analyses gener- 
ally provide good descriptions for multi-phase cooling flows 
in the cores of clusters (Thomas, Fabian, & Nulsen 1987; 
Allen, Fabian, & Kneib 1996). Since no significant radio 
emission has been detected from NGC 3923 (e.g. Birkinshaw 
& Davies 1985) it is unlikely that the hot gas is supported 
by magnetic fields. It should also be added that the p ~ r~ 2 
model appears to fit the radial profile better than the Hern- 



quist model for r ~ 100" — 300", but the effect is not highly 
significant because of the low S/N in this regio nQ 

If we fit the M oc Lr model and add dark matter fol- 
lowing the p ~ r~ 2 model, the fits are marginally improved 
over the single-component case with a minimum \ 2 = 13-3 
obtained for a ratio Mdm / Mr = 10 and Mum /Mr > 3 
(90% confidence). The ellipticity of the dark matter, 6dm, 
is essentially that of tmass for Mdm /Mr = 10 but increases 
systematically so that 6dm ~ tmass +0.10 at the lower limit 
Mom /Mr — 3. Here it should be understood that Mdm in- 
cludes only dark matter that is distributed differently from 
Lr and that the \ 2 values from the radial profile fits for 
given ratios Mdm /Mr depend on the assumed temperature 
profile. 

The Hernquist profiles are only affected when 
Mdm /Mr ~ 1 at which point the fits are improved to min- 
imum x 2 ~ 21 with large a s ~ 400". The dark matter in 
these models is required to be extremely flattened. 

Similar to our analysis of the j3 model (§@), we find 
that adding a discrete component with X-ray emission pro- 
portional to Lr does not improve the fits for the p ~ r~ 2 
model; we place a slightly stricter limit of fh g / fdisc > 4.5 



^1 We remove a larger region for the PSPC because of its larger I'll Note that the limits on e ma ss are essentially unaffected when 

PSF. the central bins are excluded. 
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(90% confidence). The Hernquist models are actually im- 
proved by the addition of a discrete component such that 
Xrod — 1-2 for fhg/fdisc = 4.6. However, essentially all of 
this improvement takes place for the central bins discussed 
above; i.e. if we exclude those bins from the fits, then the 
discrete model does not improve the fits significantly. Hence, 
a discrete component distributed like Lr does not improve 
fits to Yi x with the p ~ r~ 2 model, and the improvement 
observed for the Hernquist model occurs in the central bins 
for which other physics (e.g, multi-phase cooling flows, mag- 
netic fields) may equally affect the models. 

In Table [?] we give the spherically averaged masses cor- 
responding to the models in Table []. We only quote masses 
at a few interesting radii to emphasize that without precise 
constraints on the temperature profile (and only an aggre- 
gate constraint on the isophote shapes) we do not really con- 
strain the detailed mass profile. (Note that the spectral con- 
straints in §^ which indicate an approximate isothermal gas 
apply only for r 15/ifQ 1 kpc.) The masses of the p ~ r~ 2 
and Hernquist models agree very well within Wh^ kpc as 
expected because over much of this region the large scale 
lengths of the Hernquist model imply an approximate loga- 
rithmic slope of -2. Assuming the gas to be isothermal out 
to r > 50/ifQ 1 kpc, the p ~ r~ 2 model, with its flatter slope, 
has approximately twice the mass of the Hernquist model at 
that distance. The gas mass is less than 1% of the gravitat- 
ing mass over the entire radius range investigated. 

Using the total B-band luminosity from Faber et al. 



2 x 10 43 erg cm -2 



( M89J ) scaled to 30 Mpc, L B = 3.45/i 70 
s~ = 6.95/if 2 x lO lo Z/0, we have that the B-band mass- 
to-light ratio in solar units is ~ 5/i7oMq/Lq at 10/i^q 1 kpc 
for the p ~ r~ 2 and Hernquist models. At r = 50/i^q 1 kpc, 
the p ~ 



model has M/Lb ~ 32h,7o Mq/Lq while the density p ■ 



"Geometric Test" result is conservative since it only consid- 
ers signatures of dark matter that are distributed differently 
from the optical light. 

Although the evidence for dark matter from the Geo- 
metric Test is marginal, the results from models which em- 
ploy an explicit solution of the hydrostatic equation assum- 
ing an isothermal gas (which is supported by the PSPC spec- 
trum - indicate that M tx L is highly inconsistent with 
the radial profiles of the PSPC and HRI data (x? e d — 3.5 
for 16 dof). This particular discrepancy arises because L is 
too centrally concentrated: the derived scale length of the 
gravitating matter is approximately 1.5-2 times that of L. 
The ellipticities predicted by this M tx L model fall be- 
low the PSPC data at a significance slightly greater than 
the 90% level. Modeling the gravitating mass with a density 
run p ~ r~ 2 or with a Hernquist profile we find that the 
ellipticity of the gravitating matter is, e ma ss — 0.35 — 0.65 
(90% confidence), which is larger than the intensity weighted 
optical ellipticity (e) = 0.30. 

This evidence for dark matter which is more flattened 
and more extended than L is similar to our conclusions from 
previous X-ray studies of two other ellipticals, NGC 720 and 
NGC 1332, but at somewhat smaller significance level than 
for NGC 720 (e.g. Buote & Canizares 1997b). These results 
are consistent with analyses of known gravitational lenses 
(e.g. Keeton, Kochanek, & Falco 1997), two polar ring galax- 
ies (Sackett et al. 1994; Sackett & Pogge 1995), and flaring 
disks in spiral galaxies (e.g, Oiling 1996). The ellipticities 
of the gravitating matter derived from our X-ray analyses 
and these other methods are consistent with those of halos 
produced by CDM simulations (e.g. Dubinski 1994). 

If an isothermal gas is assumed then models with matter 



are favored over Hernquist models (and sim- 



Hernquist model has M/L B ~ l7h 70 M Q /LQ. These M/L B il ar models like the universal CDM profile of Navarro et al. 



values within 10h^ kpc are very consistent with stellar dy- 
namical studies of ellipticals (e.g. van der Marel 1991). The 
larger (extrapolated) M/Lb values at 50/ifg 1 kpc agree with 
previous X-ray studies of ellipticals (e.g. Sarazin 1997) indi- 
cating that NGC 3923 has a mass distribution typical for a 
galaxy of its luminosity. 



9 CONCLUSIONS 

We have analyzed the gravitating matter distribution of the 
E4 galaxy NGC 3923 using archival X-ray data from the 
ROSAT PSPC and HRI. Analysis of the PSPC data, which 
allows more precise constraints than the HRI data, demon- 
strates that the X-ray isophotes are significantly elongated 
with ellipticity e x = 0.15(0.09 - 0.21) (90% confidence) for 
semi-major axis a ~ 10h^ kpc and have position angles 
aligned with the optical isophotes within the estimated un- 
certainties. A bright point source located ~ 100" along the 
major axis inhibits reliable ellipticity constraints for larger 
radii. 



1997). For r ~ 100"-300" the p ~ r~ 2 model marginally fits 
the data better than the Hernquist model. However, most 
of the difference in these models occurs in the central radial 
bins where the effects of multi-phase cooling flows, magnetic 
fields, and discrete sources could affect the surface brightness 
profiles, though we have argued the effects are unlikely to 
be important (see §^). (The derived shape of the gravitating 
mass is mostly robust to these issues - Buote & Canizares 
1997b.) This support for nearly r~ 2 profiles agrees with pre- 
vious studies of gravitational lenses (e.g. Maoz & Rix 1993; 
Kochanek 1995), although a recent paper finds that density 
profiles wi th changing slo pes (e.g. Hernquist and NFW) are 
preferred ( |Williams 1997b - 

An emission component that is proportional to L can- 
not contribute significantly to the ROSAT X-ray emission of 
NGC 3923, and thus discrete sources should not affect our 
constraints on the gravitating matter (Buote & Canizares 
1997a). However, the ASCA spectral data when fitted with 
two thermal components yield a cold component, Tc = 0.55 
keV, and a hot component, Th ~ 4 keV, where th e relative 
flux of cold-to-hot is ~ 1.9 in the ROSAT band puote & 



By applying a "Geometric Test" for dark matter, which 



essentially compares the shapes of the observed X-ray 
isophotes to those predicted if mass traces the optical light 
L (independent of the poorly constrained temperature pro- 
file of the gas), we found that the ellipticity of the PSPC 
X-ray surface brightness exceeds that predicted by the con- 
stant M/L hypothesis at the 80%-85% confidence level. The 



Fabian 1997). The conventional interpretation of the hot 
component (e.g. Matsumoto et al. 1997; Loewenstein & 
Mushotzky 1997) is that it arises from discrete sources. But 
our analysis (§|5|) shows that ~ 35% of the 0.5-2 keV emission 
cannot be distributed like the optical light which would be 
expected of discrete sources. Hence, either the emission from 
discrete sources is not distributed like L, or the hot compo- 
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Table 7. Gravitating Mass and Gas Mass 

Model ^^70 ^^70 50/i^q 1 kpc 

p ~ r" 2 (1.1 - 2.8) x 10 11 /if 1 Af Q (2.4 - 6.0) x lO 11 ftf o 1 M (12.2 - 30.5) x 10 11 hf 1 M Q 

Hernquist (1.3 - 3.3) X 10 11 M Q (2.6 - 7.1) X IO^H^Mq (6.0 - 18.1) X 10 n hj^M Q 

Gas 0.14 x lO 9 /i~ 5/2 M 0.45 x 1O 9 /i~ 5/2 M 6.4 x 1O 9 /i~ o 5/2 M 

90% confidence values of the gravitating mass corresponding to the oblate and prolate models in Table 
which include 90% uncertainties in the (isothermal) temperature. The statistical errors on the gas mass 
are less than 10%. 



nent obtained from the spectral fits cannot be entir ely du e 
to discrete sources as suggested by Buote & Fabian (1997). 

The constraints we have obtained for NGC 720, NGC 
1332, and now NGC 3923 from analyses of their X-ray 
isophote shapes and radial surface brightness profiles pro- 
vide an initial demonstration of the power of X-ray analysis 
for probing the shape and radial distribution of gravitating 
matter in early-type galaxies. The next generation of X-ray 
satellites, particularly AXAF and XMM, have the capability 
to accurately map X-ray isophote shapes and orientations 
from the cores (r ~ 1") out to 10s of kpc for many galax- 
iea I The spatially resolved spectra provided by these fu- 
ture missions will allow more precise constraints on temper- 
ature gradients and the contribution from discrete sources. 
Thus, unlike most other methods, obtaining interesting X- 
ray constraints on the shape and radial density profile of 
the gravitating matter will be possible for a large sample 
of early-type galaxies since the X-ray analysis is applicable 
to any isolated early-type galaxy whose soft X-ray emission 
(~ 0.5 — 2 keV) is dominated by hot gas. 
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